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Abstract 

The Cumulative Drag Index defined recently by Prasanna [1] has been gener¬ 
alised to include the centrifugal acceleration. We have studied the behaviour of the 
drag index for the Kerr metric and the Neugebauer-Meinel metric representing a 
self-gravitating rotating disk and their Newtonian approximations. The similarity 
of the behaviour of the index for a given set of parameters both in the full and 
approximated forms, suggests that the index characterises an intrinsic property of 
spacetime with rotation. Analysing the index for a given set of parameters shows 
possible constraints on them. 
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1 Introduction 


The phenomenon of rotation plays a very important role in almost all classes 
of objects that encompass our physical universe. Particularly in the discussion 
of Inertia, rotational features characterise global effects on local physics as im¬ 
plied by Mach’s principle. Recently, Prasanna [1] has dehned a new parameter 
called the Cumulative Drag Index for stationary axisymmetric spacetimes, us¬ 
ing the notion of inertial forces within the framework of general relativity. The 
index, dehned for particles in circular orbit along the trajectory on which the 
centrifugal acceleration is zero, characterises the intrinsic feature of rotation 
through the drag induced on both co-rotating and counter-rotating particles. 
However, for practical applications, it would be useful to generalise the drag 
index to include the centrifugal acceleration also. A few years ago, astronomers 
discovered two co-spatial stellar disks in the galaxy NGC-4550, one orbiting 
prograde and the other retrograde with respect to the galactic nucleus, in the 
core of the Virgo cluster [2,3]. Bicak and Ledvinka [4] have tried to construct 
sources for the Kerr geometry using counter-rotating thin disks. If one consid¬ 
ers the galactic nucleus as a black hole, one can then use the Kerr geometry 
for the outside and have counter-streaming jets outside the ergo-region. The 
presence of such co- and counter-rotating particle streams may perhaps be 
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characterised through the drag index, dehned as 


^ ^ (J^cf + ^co - -^gr) 

(•^cf + -^co + -^gr) 

where J?cf, and jFgj,, denote, respectively, the centrifugal, the Coriolis 
and the gravitational accelerations acting on a particle in circular orbit in a 
stationary, axisymmetric gravitational field. Within the framework of general 
relativity this dehnition is unique, when one considers the spacetime expressed 
in the conformal 3+1 splitting with the four acceleration a* being expressible 
covariantly as [5,1] 

a, = + r'^VkUi) + . (2) 

The various quantities on the r.h.s. of eq. (2) are as described below: n* is the 
vector field corresponding to the zero angular momentum observers expressed 
in terms of the Killing vectors ? 7 * (timelike) and (spacelike) as 

n* = e'^(?7* + ), uj =-{v,0, (3) 

and (j) is the scalar potential 

0 = -^ H-{V: v) - v) - 0 ) • ( 4 ) 

r* is the unit spacelike vector orthogonal to n* along the circle depicting the 
orbit of the particle with a constant speed V, and 7 (= 1/Vl - ^ 2 ) is the 
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Lorentz factor. The particle four velocity, f/*, is thus expressible as 

= 7(n* + Vt^) , 

and is also equal to A(r/^ + with A the redshift factor dehned as 

= -((?], 7]) + 2n{^, r]) + , (5) 

n being the angular velocity, hlr* = e'^(r2 — f* = is the vector de¬ 

hned on the conformally projected 3-space having the positive dehnite metric 
hik = dik + niUk and V* is the covariant derivative with respect to hik = e^'^hik- 

As shown earlier, for the metric 

ds^ = {gudt^ + 2gt^dtd(t) + g^<f,d(j)^) + {grrdr^ + geedO^) ( 6 ) 


the accelerations are given as 


r..2 


Gravitational: i^gr)i = < In 


9t(f> dttgcjxj) 
94>i> 


(7) 


Coriolis: (J^co)i = 'y'^Vn^iVjTi - ViTj) 

A ^)g<p<t)di{^gt(p/g,f)(j >), 


( 8 ) 


Centrifugal: (^cf)i = ( 7 ^^) 7*Vfcf, 

A‘^{Q-u) 


-g^^di < In 




9t<f) 9tt9<p<f> 


(9) 
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2 Kerr spacetime 


Taking now the specific example of Kerr spacetime 

\ , , 4mra . 9 ^,,, 5. 9,,9 

] dt -—— sm 6dtd(j)+ — sm 

J Zj 

( 10 ) 

where B = (r"^ + a?Y — Aa^ sin^ A = — 2mr and S = cos^ 9, 

and considering a particle in circular orbit on the equatorial plane {6 = 'k/2), 
it can be seen that the index is 

C = |m[r^(r — 2m) + 2a^r{r^ — 8mr + lOm^) + a^(r — 6 m)] 

+2amf2[r‘^(r + 2m) + 2a^r(r^ + 4mr — lOm^) + a‘^{r + 6 m)] 

+f2^(r^ + a^r + 2ma^)[r^(r — 4m) + a^r{r‘^ — 5mr + lOm^) — 3a"^m]| 

j |A(r^ + a^r + 2 a^m)(—m(l — afl)^ + f 2 ^r^)|. ( 11 ) 

It is clear that, of the two inhnities of the index, one appears at the event 
horizon (A = 0) while the other depends on both a and f2 and appears for a 
given a and f2 at r = [m(l — Fig. (1) shows the nature of C at the 

three locations, rp^p (the prograde photon orbit), rphr (the retrograde photon 
orbit) and rcfo (the orbit where the centrifugal force is zero). 

While at rphp the index is positive only for a very small range of for counter¬ 
rotating particles (f 2 < 0 ), at rphr the index is positive for the same range of 
f2 for co-rotating particles only. On the other hand, as was discussed earlier 



5 



in [1], at Tcfo the index is positive for both co- and counter-rotating particles, 
but again for a very narrow range of values of hi (Fig. (Ic)). This change of 
behaviour of C at the two photon orbits arises due to the following reason: 
When the centrifugal force is not zero, the two zeros of the denominator of C 
outside the event horizon corresponding to the fixed value of r and a, are at 


ffi = 


m 


aWm + r 


3/2 


and 


= 


a\/m 


. 3/2 


( 12 ) 


and corresponding to these two hi values the numerator of C factors as 

2mr^^'^{—2ay/m + — r^'^^)[4a^mr — (a^ -|- r^)^] (13) 


and 

2mr^^‘^{2ay/m + — r^^‘^)[Aa‘^mr — (a^ + r^)^], 

respectively. Thus the zero at hli cancels with the numerator at the prograde 
photon orbit, while the one at 112 cancels with the numerator at the retrograde 
photon orbit. 

For many applications, often one takes the view that the linearized Kerr metric 
might be sufficient to incorporate the relativistic effects, when the body is 
slowly rotating. In order to examine this, let us consider the nature of C under 
this approximation. The three accelerations acting on a particle in circular 
orbit, approximated to terms linear in the Kerr parameter a, are given by 
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T ^ 

"■■“rH r) ’ 


T = 

'j CCi 


6amQ 


(r^ — 2mr — ’ 


^cf = 


r2(r — 3m)(f2^r® — f2r^ + 2mQr^ + 4mar — Sam^ 
r(r — 2m){r — 2m — Vfir'^Y 


( 14 ) 


and thus the index is 


IC = 


-m{r — 2m) — 2amVLr — VL^m{r — 5m) 


+6amf2^r^ + — 4m) /{r — 2m) 


m{r — 2m) — 2amVLr — f2^r^(r — m) + 6amQ^r^ + 


(15) 


The very hrst change one notices is that, neglecting and higher order terms 
in a, moves the inhnity at the horizon to r = 2m, as this would now represent 
the horizon, just like in the static case. Similarly, the orbit where the centrifu¬ 
gal acceleration is zero also coincides with that in the static case, viz., r = 3m, 
for all f2. Fig. (2) shows the index for the linearised version at the two photon 
orbits and at the orbit on which the centrifugal acceleration is zero. Com¬ 
parison of these plots with those for the unapproximated C (Fig. (1)) clearly 
shows that for particles with angular velocity |f2| > 0.3, the behaviour is ex¬ 
actly same with or without the approximation at the photon orbits whereas 
at the orbit with J^cf = 0, the similarity is striking for all values of 


As it would be almost impossible for particles to have a low value of close 
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to photon orbits (as they would be relativistic), the behaviour of the index 
shows that the linearisation approximation for the forces is amply justihed 
for all practical purposes. However, as a matter of principle one hnds that for 
very low values of hi the behaviour of the linearised version differs from that 
for the exact version, the difference arising mainly because of the centrifugal 
acceleration being non-zero. If one considers the behaviour of ^C as a function 
of r for hxed a and H, it seems to be exactly like for C, the expression without 
approximation. 

Looking at the overall feature of the inertial accelerations it then seems that 
for understanding the frame dragging coming from rotation, for practical pur¬ 
poses of considering the forces, it may indeed be sufficient to calculate the 
gravitational, Coriolis and centrifugal acceleration in the linearised approxi¬ 
mation, as given in eq. (14). 

It is further interesting to consider the Newtonian limit of the accelerations 
with lowest order corrections in the centrifugal acceleration, as given by: 



With these expressions the drag index turns out to be 

(r^ — 3mr^)r2^ -|- 2amVL + m 

NC = - - - -. 

(r^ — 3mr‘^)Q‘^ + 2amQ — m 
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( 16 ) 


(17) 



Fig. (3) shows the plots of NC as a function of ^2, for hxed r and a (3a,b) and 
as a function of r for hxed a and f2 (3c,d). As the Newtonian approximation 
can be valid only for larger values of r, it is clear that the index is positive 
for both CO- and counter-rotating particles for |f2| > 0.1, independent of the 
values of a. 

3 Disk spacetime 


We shall next consider the behaviour of the index for the Neugebauer-Meinel 
metric [ 6 ] for a rigidly rotating disk of dust as given in the discussion of its 
dragging effects by Meinel and Kleinwachter [7]. The metric components and 
their hrst radial derivatives of signihcance are 



where is the angular velocity of the disk, 2 ; represents the relative redshift 
of photons from the centre of the disk measured at inhnity, /i is a parameter 
dehned through the relation p = 2 f 2 ^rQ(l -|- z)'^ and tq represents the rim of 
the disk. Using eqs. 1, 7-9 and 18, the drag index may be obtained as: 


c = + 2/i{4(/i - 1 ) 12 ^ - 5/if2dfi + 

+2{(4 - 15/i + + 3p(4 - 5p)fldf2 + /i(6p - l)n‘^}z^ 

2{2(4 - 9p + 4p2)f2^ - (8 - 26/i + lh^i^)VLdVL + p(7/i - 8)f22};23 
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+{Vt-Vtdf{S-U^i + b^j?)z^ 


fiiyid — ^)(l “I" 2 :){r 2 (i(l -\- z) -\- f2(l — ^)}{/i(l + z')^ — 2^:^} , 


(19) 


where is the angular velocity of the particle. 


It is clear from the above expression that the two zeros of the denominator, 
where the index blows up, correspond to the two circular geodesic orbits, 
as the sum of the forces acting on the particle is zero for these parameters. 
While the zero at 11 = H^, corresponds to the prograde geodesic, the one at 
H = —lld(l + z)/{l — z), corresponds to the retrograde geodesic as also shown 
by Meinel and Kleinwachter. However, one hnds that if 2 : > 1, the second zero 
occurs at a positive value of H as shown in the Fig. (4d). Fig. (4) shows the 
plots of the index for four different values of /x. 


4 Discussion 


The presence of centrifugal acceleration does bring in a difference in the be¬ 
haviour of the drag index at the two photon orbits, with the co-rotating ones 
having a positive value for a narrow range of H at the retrograde photon orbit 
and the counter-rotating ones having a similar feature at the prograde photon 
orbit. However, if the black hole is slowly rotating (a 1, negligible), then 
adopting the linearised version of the acceleration changes the behaviour of the 
index only for very low values of H (|11| < 0.3) at the photon orbits, whereas 
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for higher values of hi the behaviour resembles that of the full C without any 


approximation. On the other hand, for given a and O, as a function of r the 
radial distance parameter, the index shows no change with approximation, 
thus indicating that the drag index signifies something intrinsic to the space- 
time with rotation, as its behaviour for both co- and counter-rotating particles 
appears similar, from the point of view of a locally non-rotating observer. 

In order to appreciate the significance of the index defined, one can consider 
its behaviour for the physical case of a rotating disk as depicted in Fig. (4). 
Comparing with the discussions of Meiner and Kleinwachter, we find that 
the constraints on the value of appearing through the parameter p, is well 
reflected in Figs. (4a,b,c). While Fig. (4a) shows clearly the existence of two 
geodesic orbits (pro and retro) for p = 0.1, Fig. (4b) shows for p = 0.5, 
the existence of only a prograde geodesic. In fact it is instructive to compare 
this with Fig. (lb) which corresponds to C at the retrograde photon orbit for 
Kerr, and thus conclude that p = 0.5 corresponds indeed to the last possible 
retrograde geodesic orbit. Further, from the zeros of the denominator of C 
(eq. 19), the appearance of the retrograde geodesic corresponds to the particle 
angular velocity = —fld(l + z)/{I — z). In fact, for this value of the 
numerator of the index factors as 


8 (2/i - l)nlz{l + zf{^^-z- ^^z)/{l - zf , (20) 
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which clearly shows that the zero of the denominator cancels with the term 
{2/j, — 1) in the numerator for jj, = 1/2, exactly similar to what happens in 
the case of the Kerr metric. It is obvious that for the case z > 1, i.e., fi > 
1.3519, this value of becomes positive, the infinity of the index appearing 
for a prograde orbit. However, this orbit turns out to be spacelike {v > c). 
Thus, it appears that the counter-streaming particles would have a limitation 
from the point of view of their redshift and this could play an important 
role in the analysis of particle motion in the disk associated with the Virgo 
cluster. From this discussion it follows that given any stationary axisymmetric 
metric representing an astrophysical situation, one can straightaway determine 
constraints on possible physical parameters characterising the geodesic orbits 
through the behaviour of the drag index. 

Though Newtonian physics does not directly predict anything regarding the 
nature of spacetime as influenced by rotation, it is amazing to see that the cu¬ 
mulative drag index shows exactly similar behaviour in the Newtonian approx¬ 
imation for the case of the full Kerr geometry as well as for the Neugebauer- 
Meinel disk geometry (Fig. (5)), for all values of a and Qd, either prograde or 
retrograde. Thus it is clear that the cumulative drag index defined above char¬ 
acterises an intrinsic property of spacetime with rotation, which goes beyond 
approximations. The fact that it is positive for both co- and counter-rotating 
particles having reasonable angular velocities, outside the ergo-region of a 
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black hole or in a self-gravitating disk, clearly supports the possibility of sus¬ 


taining counter-rotating streams. In fact, this analysis points out a constraint 
on the value for counterrotating streams to be < 0.285, which might be tested 
in the case of streams encountered in the Virgo cluster. Eventhough the drag 
index itself does not measure any observable quantity directly, it characterises 
intrinsic rotation for stationary, axisymmetric spacetimes, be it in an empty 
region (like Kerr geometry) or within a rotating disk (Neugebauer-Meinel class 
of solutions), independent of approximations. If the metric potentials depend 
upon a directly observable parameter then the behaviour of the index for dif¬ 
ferent combinations of rotational parameters could possibly give constraints 
on the physical parameter, which may be measurable. Further, as the free 
orbits are dehned through the equilibration of the forces acting on a particle, 
the index would go to inhnity and thus studying the behaviour of the index 
in general yields information for orbits both free and otherwise. Thus for the 
case of stationary axisymmetric metrics, the location of geodesics for given 
angular velocity and rotational parameters may be identihed by plotting the 
drag index without necessarily solving the equations of motion. 
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Figure Captions 


Figure 1: 


Figure 2: 


Figure 3: 


Figure 4: 


Figure 5: 


C (cdi) as a function of (a = 0.5m) at the prograde photon orbit (a), the 

retrograde photon orbit (b), and the orbit with Ta = 0 (c). 

iC as a function of (a = 0.5m) at the prograde photon orbit (a), the 

retrograde photon orbit (b), and the orbit with = 0 (c). 

iVC as a function of for hxed r and a (a, b), and as a function of r for 

hxed a and hi (c, d). 

C as a function of for a rotating disk of dust with (a) fi = 0.2, = 0.1, 

(b) = 0.5, = 0.1, (c) /i = 2, = 0.99 and (d) /i = 1.9, = 0.1. 

Comparison of C and NC for Kerr (a, b) and the rotating disk (c, d). 
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at rphp (a), rphr (b) and r^o (c). 


a function of (a = 0.5m) 
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Fig. 5. Comparison of C and NC for Kerr (a, b) and the rotating disk (c, d). 
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